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Abstract 

We present a method that allows us to find the asymptotic form of various char- 
acteristics of disordered systems in the strong localization regime, i.e., when either 
the random potential is big or the energy is close to a spectral edge. The method 
is based on the hypothesis that the relevant realizations of the random potential in 
the strong localization regime have the form of a collection of deep random wells 
that are uniformly and chaotically distributed in space with a sufficiently small den- 
sity. Assuming this and using the density expansion, we show first that the density 
of wells coincides in leading order with the density of states. Thus the density of 
states is in fact the small parameter of the theory in the strong localization regime. 
Then we derive the Mott formula for the low frequency conductivity and the asymp- 
totic formulas for certain two-point correlators when the difference of the respective 
energies is small. 

PACS numbers: 05.60.Gg, 72.15.Rn, 72.80.Ng 



1 Introduction 

It is widely accepted and proved rigorously in many cases that elementary excitations in 
disordered media are localized if the disorder is strong enough or/and the energy of the 
excitations is close enough to the band edges. The idea dates back to the famous paper 
0] by P. Anderson who emphasized, in particular, the aspects related to the transition 
from localized to delocalized states. The idea was further developed by N. Mott and I. 
Lifshitz (see e.g. their review works [213 ED])- In particular, it was I. Lifshitz who singled 
out the regime of high disorder or low energy where the localization is most pronounced. 
This regime is now known as the strong localization regime. According to I. Lifshitz, in 
this case, the pertinent realizations of the random potential have the form of a collection 
of deep potential wells which are so rare and whose form is so irregular that the quantum 
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mechanical probability for tunnelling through a macroscopic number of the localization 
wells vanishes. 

The study of localization and relevant physical characteristics of disordered systems 
can be reduced to the study of moments of the density operator p E = S(E — H), where 
H is the (one-body) Hamiltonian of the system. By using the coordinate representation, 
we can write the /-th moment (7-th correlation function) as follows: 



where the (...) denotes averaging with respect to the disorder. 

The simplest case of the correlation function corresponding to I = 1, x% — y\ = x: 



i.e., to the average of the local density of states Pe(x, x), is known as the density of states 
(DOS) of the system. 

I. Lifshitz suggested a non-perturbative method of computing the asymptotic form 
of the DOS in the strong localization regime JH]- The above description of typical re- 
alizations of the random potential is implemented in this method by the assumption of 
independent quantization of a quantum particle in each localization well (see [2TH |2~H 17]). 
thus the complete localization of a particle in an exponential neighborhood of each well. 
A rigorous proof of the complete and exponential localization in the strong localization 
regime was given by J. Frohlich and T. Spencer |§1 12*2*1 12*7]. 

In both these important results of the localization theory the tunnelling between the 
localization wells plays no significant role. In Lifshitz's argument, other wells are sim- 
ply ignored. The crucial ingredient in the rigorous proof of the complete localization 
in the strong localization regime is a rather sophisticated probabilistic extension of the 
Kolmogorov- Arnold- Moser theory (known as the multi-scale analysis) which allows one to 
verify that tunnelling between wells is strongly suppressed, and therefore does not change 
qualitatively the picture, suggested by the independent wells quantization assumption. 

The DOS determines equilibrium properties of a disordered system in the one-body 
approximation, i.e., of the ideal Fermi gas in a random external field. The study of kinetic 
properties of the gas and of interaction effects requires the knowledge of higher moments 
of the density operator p E , especially the second moment K%. Important quantities 
that can be expressed via K 2 are the density-density correlator and the current- current 
correlator ^UJ|^]. These correlators allow us to answer relevant questions concerning the 
nature of localization and the behavior of the conductivity and other physical character- 
istics. 

The complete localization of states in a certain interval of energies implies that the 
zero temperature d.c. conductivity vanishes if the Fermi energy lies in this interval (see [2] 
for a proof and a discussion). On the other hand, since the energies of localized states are 
dense, the zero-temperature a.c. conductivity is expected to be non-zero for any non-zero 
frequency v of the external field. It was N. Mott who first proposed "resonant" tunnelling 
between pairs of wells as a mechanism of the low frequency a.c. conductivity in localized 
systems j23j. According to Mott, one can view those states, resulting from independent 
quantization in each localization well (localization center in Mott's terminology), as a 
kind of "bare" states. They decay exponentially in the distance from the corresponding 
localization center. Two (several) bare states with widely spaced centers but with suf- 
ficiently close energies can "resonate". This leads to the two-center states (resp. multi 



Ki(xi, ...,x t ; y h ...,yr, E x , E t ) = (p El (xi, yi)-.p El {x h yi)), 
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center states), whose energies are exponentially close in the separation between the cen- 
ters. The condition for a pair of wells to be in resonance determines the distance between 
resonating wells, thereby determining the characteristic value of the dipole moment of two 
bare states of wells, and the square of the dipole moment is, in essence, the conductivity 
according to the linear response theory (see formula (|3.4jl below). This observation leads 
to the following asymptotic expression for the low frequency conductivity: 

a(u,E F ) = Ap\E F y (log^) d+1 (1.3) 

in the case, where 

T<z/<£ F . (1.4) 

Here T is the temperature, v the frequency of an alternating external field, E F the Fermi 
energy (supposed to be in the localized spectrum). A and u are determined by the 
fundamental constants and by the random potential. 

Formula (jl.3j) was discussed in many works (see e.g. [7J |Hl 13 EH HU HTJ 123 G2] and 
Section 5). However, a consistent "first principle" derivation of the formula is still not 
available in a general multi-dimensional case. We mean a derivation based on the Kubo 
formula (see formulas (J3.2j) - (J3.4j) below), in which the two-point correlation function is 
computed for a given random potential in the asymptotic regime (11.4)1 . 

The fact that such a derivation is still missing encourages us to present in this paper 
a heuristic method that allows us to obtain formula (jl.3)) and some other two-point 
correlation functions (i.e., (|l.l)l for I = 2 and \Ei — E 2 \ := v <C \Ei\, \E 2 \), and that, 
we believe, clarifies Mott's initial arguments. 

The method is based on the above hypothesis on the form of pertinent realizations of 
the random potential as systems of deep and rare localization wells. Viewing the density 
of wells as a small parameter of the theory, we apply a version of the virial expansion to 
compute the leading contribution to the moments Ki of (jl.lj) for 1 — 1,2. In particular, 
by applying this procedure to the DOS, we find that its leading order is the density of the 
localization wells. This shows that the small parameter of the theory is the DOS itself, 
whose smallness is known to be an important condition for localization. Furthermore, 
we find that the leading order of the pair correlation functions, the a.c. conductivity in 
particular, is determined by two-center states, resulting from resonant tunnelling between 
a pair of localization wells, in agreement with Mott's ideas. This leads to formula ()1.3|) 
and, therefore, supports the idea of pair approximation in Mott's derivation of 1)1.3)1 . 
Among our other results, we mention high peaks of some pair correlation functions (see 
(14.1)) and ()4.2)1 below), appearing in a neighborhood of the origin and on the "resonating" 
distance, determined by the frequency of the external field. Analogous peaks were found 
before in the one-dimensional case for strong localization as well as in the weak 
localization regime ^T]. However, in these cases, the peaks are of the order p 2 (E F ), while 
in the general (i-dimensional case, the peaks are of the order p 2 (E F ) (\ogv G /v) ~ 1 , i.e., 
much bigger in the regime (jl.4j) (see also ^| for a similar result). 

The paper is organized as follows. In Section 2 we outline the method. In Section 3 
the Mott formula ()1.3)1 is derived. In Section 4 we derive asymptotic formulas for binary 
correlators and in Section 5 we comment on our results and on their relations to known 
results. 
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2 Method 



2.1 Effective potential 

It was already mentioned in the introduction that extensive studies of the strong local- 
ization regime show that the phenomenon is determined by realizations of the random 
potential, containing deep and rare potential wells. For a potential unbounded below (like 
the Poisson potential (|2.6|) below) the large parameter of the theory is the absolute value 
of the energy and/or the amplitude of the potential. These two cases of the strong local- 
ization regime are manifestations of the simplest mechanism of localization: capturing a 
quantum particle in strong and rare fluctuations of a random potential 1 . 

In other words, for an overwhelming majority of eigenf unctions ipj, corresponding to 
the strongly localized part of the spectrum, there exists a point the center of the 
localization well, such that r/> - decays as exp{ — \x — Here r$ is the localization 

radius of ipj. The localization centers have to be uniformly and chaotically distributed in 
space and the distances between them have to be much bigger than the typical localization 
radii and than the radii of the localization wells. Hence, one has to expect an effective 
"decoupling" between the localization wells. 

One obtains a simple form of this picture of the strong localization regime by replacing 
the random Schrodinger operator by the direct sum of operators, each of them defined 
in a certain cell, containing a single localization well. This procedure of independent 
quantization in isolated cells is supported by and even instrumental in studies of the den- 
sity of states, the interband light absorption coefficient, and other spectral and physical 
characteristics of disordered systems (see e.g. [3 [HI 121]), as we ^ as °f the probability 
distribution of spacings between adjacent energy levels (see (211123) in the strong localiza- 
tion regime. However, the procedure is not appropriate in studies of transport properties 
of disordered systems. This is why we replace the procedure of independent quantization 
in isolated cells by the less restrictive assumption, according to which relevant properties 
of the strong localization regime can be described, assuming that any shortly correlated 
and smoothly distributed random potential can be replaced by an (effective) potential of 
the form: 

3 

Here {£ •} are the Poisson random points of the density fi, modelling the centers of the 
localization wells, and the random functions {vj} are independent of each others and 
independent of the {£_,•}. The {vj} model the shape of the localization wells. We assume 
that all v/s have a finite range and the typical radius a of v/s is related to the typical 
distance between wells as 

a < ^ 1 ' d . (2.2) 

The density /x of the localization centers is not known and has to be found self-consistently. 
The density as well as the shapes of the wells may depend on the energy interval in 
question. 

In other words, we believe that the strong localization regime possesses a certain 
robustness (insensitivity) with respect to a concrete form of random potential, provided 

■"■We mention another localization mechanism : enhanced backscattering. The mechanism is respon- 
sible for localization at high energies in the one-dimensional case, and for weak localization effects in 
arbitrary dimensions. 
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that it is translation invariant in the mean, shortly correlated, and smoothly distributed 
(the last two properties facilitate the localization because they make it more unlikely that 
different localization wells are of the same shape, thereby suppressing tunnelling between 
different localization wells). One may say that our ansatz (|2.1|) replaces impenetrable 
walls between cells of the independent quantization procedure by a kind of "soft" walls, 
that strongly suppress particle mobility but do not exclude it completely. 

To avoid technicalities, we will choose a simple form of the localization wells, setting 

Vj(x) = g 3 v (y/gjx) , (2.3) 

where v(x) is a finite range potential well and {gj} are independent identically distributed 
random variables, independent of {£ -} and assuming arbitrary big positive values accord- 
ing to a smooth probability density p(g). 

Summarizing, we can write the following formula for the effective potential 

Veff(x) = Y^gjv {y/gj(x - Q) . (2.4) 

3 

It should be noted that similar random functions are widely used in localization theory 
as "bare" random potentials in the Schrodinger equation (see e.g. [Tfi| l21j). We mean the 
potentials of the form 

u(x — Xj), (2.5) 

i 

where u is a non-positive function of a finite range (the single-impurity potential). In 
the case, where {9j} are independent identically distributed random variables and {xj} 
form a regular lattice, the potential models a substitutional alloy, and in the case, where 
9j = 9 = const for all j and {xj} are completely chaotic (Poisson) random points of the 
density c, the potential 

u(x — Xj), (2.6) 

3 

models an amorphous medium. Assuming that c is large, 9 is small but c9 2 = D is fixed 
and shifting the energy by the mean value 




of the potential ()2.6|) . we obtain a Gaussian random potential with zero mean and with 
the correlation function 

D J u(x- y)u(y)dy. 

In a more general case, where the {xj} are completely chaotic and the {9j} are identically 
distributed random variables, independent of each others and of the {xj}, (|2.5|) is a 
generalized Poisson potential. 

We would like to stress here that while our effective potential f!2.4|) is similar to a 
generalized Poisson one (because of random g/s), these two should not be identified. In 
particular, the density c of the impurity centers {xj} in (12 .5|) is not the density fi of the 
localization centers {£ •} in (12. 4|) (fi is usually is much smaller than c), and the functions 
9jU in (J2.5J) . modelling the single impurity potential, have little in common with the 
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functions Vj in (J2.4j) . modelling the form of the localization wells. The latter are formed 
by sufficiently large and dense clusters of impurities in which the inter-impurity distances 
are much smaller than the typical distance c~ l / d between impurities centers {xj}. For 
example, if the "bare" random potential is given by ()2.6j) . then it can be shown than the 
number of x/s in a typical localization well is of the order logE /u(0) >> 1 [2~Tj . 

2.2 Density expansion 

Recall that an important property of the effective potential is the small density fi of the 
localization centers (cf ()2.2j0 . We describe now a technique that will allow us to use this 
property. 

Let {Fi(x\, xi)}i>q be a system of functions of / (^-dimensional variables x%,...,Xi 
(F Q is a constant). We denote the set (x\, ...,xi) as X. Suppose that the system {-Fj}j>o 
satisfies the following conditions (we do not indicate explicitly the index /). 

(i) Translation invariance: for any d- dimensional vector a 

F{X) = F(X + a), where X + a = (x x + a, xi + a). 

(ii) Additive clustering: 

F(X U (Y + a)) - [F(X) + F(Y)] -»• 0, as a -»■ oo. (2.7) 

and the decay of the l.h.s. of ()2.7|) is fast enough (it will be exponential below). 

For any system of functions, possessing these properties, we can write the combinato- 
rial identity 

w = EE(- 1 ) Wf ( z ). ( 2 - 8 ) 

YcX ZcY 

where N(X) is the number of points of X. 

We will use this identity in the case, where X, Y, Z are the sets of random Poisson 
points entering in the effective potential (|2.1|) . Recall that an infinite system {£ ■} 
of Poisson points of density /x in the d- dimensional space can be asymptotically described 
as a system of random points £ 1( ...,£ N , uniformly distributed in a cube A, provided that 
the "thermodynamic" limit N — > oo, |A| — >• oo and A/|A| — \i is carried out (we will 
denote this limiting transition by A — >• oo). By using this fact and identity ()2.8j) . we can 
write that 

hm \A\- l {F N (t 1 ,...,Z N )-F ) = fj,(Fx — F ) (2.9) 

A— >CXD 

+ ^ J[F 2 (x)-2F 1 + F }dx + ... 

where the symbol (...) in the l.h.s. denotes averaging with respect to the Poisson points 

In view of (|2.4j) . we will need a more general formula in which the role of is played 
by pairs (£,j,gj), where {gj} is a system of independent random variables of common 
density p(g) which are also independent of the The corresponding formula can be 
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obtained from (|2.9|) . written for fixed g/s and subsequently integrated with respect to 
g/s with the probability density p(g). This yields 

lim \A\- 1 (F N ((^,g 1 ),...,^ N ,g N )) - F ) = f{F 1 {g 1 )-F )fJi{g 1 )dg 1 (2.10) 

A^oo J 

+ \ J \ F 2{x] gi, g-i) - Fi(gi) - F x {g 2 ) + Fo]n{g 1 )n{g 2 )dxdgidg2 + 

where now the symbol (...) in the l.h.s. of this formula denotes averaging with respect to 
{^■} and {gj} and 

H(g) = w(g). (2.11) 
2.3 Density expansion of the DOS 

Now we apply the expansion described above to the density of states of the Schrodinger 
equation. We use the self-averaging property of the DOS, according to which 

p(E)= \im(\ A]' 1 J2^E-E n )), (2.12) 



A^oo 

n>l 



where {E n } n >i are the energy levels of the Hamiltonian H\ defined by the Schrodinger 
equation with the potential (j2.4j) in the cube A. 

Comparing the l.h.s. of ()2.10|) and the r.h.s. of ()2.12|) . we conclude that in this case 
the role of F t in (|2~TU|) play 



^25(E-ES>((x 1 ,g l ),...,(x h g l ))), 

n>l 

where {E n ((x±, gi), (xi, gi))} n >i is the negative spectrum of the /-wells Hamiltonian 



l 

H® = -A + J2 9j v iV9](x - Xj )) . (2.13) 

3=1 

Thus, applying (|2.1(J|) to the DOS and taking into account that we are interested in 
negative energies of large absolute value and that = —A has no negative spectrum, 
we find that the term p(°'(E) with / = (the zero- well contribution) is absent in the 
expansion. Hence, the leading contribution in /i to the DOS is due to the one-well term 
of the expansion: 

P (1) (£) = E / 5 ( E - Ei n ] )^9)dg. (2.14) 



For the well of the form gv (yfgx\ we have: 

E<f> = ge n , (2.15) 
where {e n } n >i are the negative eigenvalues of the dimensionless operator — A + v(x). Thus 



P 



M <*) = E"(f)irr 

n>l /ii 
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According to the spirit of our approach the density p(g) should decay sufficiently fast as 
g — > oo. Thus the leading contribution to p^(E) is due to the first term of the sum, i.e., 
we can use the approximation 

\ e iJ \ £ i\ 

Normalizing the well v by the condition 

ex = -1, (2.17) 

we can write 

p^(E)^^-E). (2.18) 

The last formula is a version of the well known "classical" asymptotic formula for the 
DOS valid for smooth random potentials. By choosing as a randomizing parameter of 
the wells Vj in (|2.1|) their ground state energies, we can show that an analogue of 1)2.16)1 
allows us to obtain also "quantum" versions of asymptotic formulas for the DOS valid for 
singular f's (see |3T] for the respective terminology and results). 

It can also be shown that the two-well contribution to the DOS is of the order 0(p?). 
We postpone the corresponding argument to Section 5.1. Thus the two- well contribution 
is negligible with respect to the r.h.s. of ()2.18j) . We conclude that the unknown (and 
small) function fi(g), determining our effective potential and having the sense of the 
probability density to find a well of amplitude lying between g and g + dg with center in 
an infinitesimal neighborhood of a given x, coincides in our approximation with the DOS 
of the Schrodinger operator. This important conclusion makes our scheme self-consistent. 
It corresponds to the basic ingredient of the Lifshitz approach, according to which the 
DOS is the probability density of the localization wells, having the ground state energy 
E |201 ■ This interpretation of the DOS is widely used in the theory of disordered systems 
[7J|^]. In our approach it is a simple consequence of the ansatz (|2.4|) and of the expansion 
formulas of the previous section. 

Let A be an interval of values of random variables gj, lying in the strong localized 
spectrum with width much smaller than typical values of the g's under consideration. 
Then Jl — f A p,(g)dg will be the density per unit volume of wells, whose amplitudes are 
in A, and pT x l d will be the typical distance between these wells. Our approach is based 
on the assumption that typical distances between wells are much larger than the typical 
radii of the localization wells (cf ()1.4j) ). In the case of the effective potential 1)2 .4]) . this 
assumption can be written as 

g- 1/2 < Ji- 1,d . (2.19) 

3 A.C. Conductivity 
3.1 Generalities 

Recall that from the point of view of statistical physics, we are dealing with an ideal gas 
of electrons in the external random field V(x) (one-body approximation). In this case, the 
linear response theory leads to the following formula for the tensor of the zero-temperature 
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a.c. conductivity of a macroscopic system of spinless electrons in an external spatially 
homogeneous electric field of the frequency v at zero temperature: 



a a p(u,E F ) = lim vre 2 |A| 1 5(E F + 

A— ioo ' J 

m^n 



E m )6(E F -E n )V^vW, 



where Vmn are the matrix elements of the velocity operator iS7 a between the states ijj m and 
ip n of the system, confined to the box A. In the case of a random potential, homogeneous 
in mean and weakly correlated, the conductivity is self-averaging [21]. Thus we have in 
the thermodynamic limit, assuming for simplicity that the system is rotational invariant 
in mean: 



a/3 



ire 

y,E F ) = —5 aP a(u } E F ) } 



a{v, E F ) = lim / |A|- X £ 6(E F + u- E m )8(E F - E n )\V mn \ 2 \ 

\ mj^n I 



where 



IK. 



a) |2 
mn I ' 



(3.i; 



(3.2) 



(3.3) 



a=l 



Since V = i[H, X], where X is the coordinate operator, we have \Vmn\ = \{E m — E n )X. 
and (j3.2J) can be written as 



mn I j 



a(is, E F ) = v 2 lim ((A)- 1 V 8(E F + u- E m )5(E F - E n )\X mn \ 2 ). 



m^n 



where 



\X„ 



Q=l 



(a) 1 2 
mn I ' 



(3.4) 



(3.5) 



Note that we keep the frequency v non-zero while making the thermodynamic limit 
A — > 00 in the above formulas. This prescription is well known in kinetic theory and 
is reminiscent of keeping non-zero magnetic field while making the thermodynamic limit 
for a ferromagnetic system in order to obtain non-zero macroscopic spontaneous magne- 
tization. Another way to obtain non-zero d.c. conductivity is to set v = in ()3.2j) but to 
replace the 5-functions by a sharp function of width 7] (usually by the Lorenzian). This 
corresponds to an imaginary shift in energies instead of a real-valued shift v (see e.g. 
j2], where the an imaginary shift is used). In this paper, we will use the formula (J3.4j) . 
assuming always that the frequency is non-zero, but small compared to the Fermi energy, 
i.e., we will assume that inequality (jl.4j) holds. 



3.2 Computation 

Now we are going to apply the density expansion formula (j2.10J) to the a.c. conductivity. 
Comparing ()3.4|) and (j2.10|) . we choose the functions Fi in this 

v 2 S(E F + v - E$)8{E F - E^IX® | 2 , (3.6) 

mj^n 
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where {En } n >i are negative levels of the /-wells Hamiltonian (j2.13j) . and 

X$ n = j x^{x)^dx, (3.7) 

.Here the {Vn^}n>i are the bound states of ()2.13|) . By the same reason as in the case of 
the DOS, the zero-well contribution to the conductivity expansion is absent. Let us 
show that the one-well contribution is also absent. Combining (jH.fi j) for I — 1 and 
(I2.15j) . we obtain 




where {E n } n >i are the bound state energies (j2.15j) of the one- well Hamiltonian H^\g) = 
—A + gv(yfgx), and Xmn is the coordinate matrix element between the corresponding 

states }. Non-zero contributions to this expression are due to the pairs (m,n) such 
that 

ge n = E F , ge m = E F + v ~ E F , g(e n - e m ) = v. (3.8) 

Denoting by e the typical value of the levels £ n 's of the potential well v and by 8e the 
typical value of the spacings |e n +i — £ n \, we see that the above conditions are incompatible 
if g5e 3> v, i.e., if Epbeje ^> v. Since dimensionless, the last condition is just 

another form of our basic condition (jl.4j) . 

The two level contribution ,Ef) to the a.c. conductivity is (cf (j2.10j) ): 

a^(u,E F ) = y E / ^ + v ~ ti$?W*F ~ E^) (3.9) 
x \ X mn\ 2 K9i)K92)dgidg 2 dy, 

where {E n } n >i are the bound state energies of the two- well Hamiltonian 

H^{{x 1 ,g l ),{x 2 ,g 2 )) = -A + g 1 v 1 +g 2 v 2 , (3.10) 

in which 

v k (x) = v (^(x -x k )), k = 1,2, 

(2) 

y = x\ — x 2 , and Xmn are the corresponding coordinate matrix elements. 

In view of our basic condition (j2.19jl . we have typically \x\ — x 2 \ 3> max g 1 l^ 2 . Hence, 
according to general principles of quantum mechanics, each level of (j3.10j) should be (ex- 
ponentially) close to a certain level of one of infinite distant wells, and each eigenfunction 
is (exponentially) close either to an eigenfunction of one of the wells (non-resonant case) 
or to a linear combination of the eigenfunctions of the both wells with coefficients of the 
same order of magnitude (resonant case). 

To make this description more quantitative, consider the one-well Hamiltonians 

Hp = -A + g k v k , k = 1,2, 

corresponding to (j3.10j) . Normalize the potential well v (x) by the same condition (j2.17j) . 
Then the lowest eigenvalues of , k = 1,2 are —g k , and the corresponding eigenfunc- 
tions are 

9k( x ) = 9k <P (V9k(x - x k )) , k = l,2, (3.11) 
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where ip(x) is the ground state of the dimensionless operator —A + v(x). The function 
ip(x) decays exponentially in x with rate 1. Hence 

I/O 

ip k (x) ~ exp (-y/gk\x - x k \) , \x-x k \^>g k . (3.12) 

Since we will be interested mostly in the resonant case, we assume that g± >2 — \Ep\, i.e., 
the radii of the ip k , k = 1, 2 in ()3.11|) are of the same order of magnitude 

g^P ~ n = \Ef\- 1 * (3.13) 

hence ri <C \x\ — x 2 \. 

In this situation we can find the lowest eigenvalues of in the framework of the 
widely used approximation, in which is replaced by its projection on the span of the 
functions (p 1 and (p 2 2 ■ The diagonal entries of this 2x2 matrix are 

{VkiH {2) ¥k) = ~9k + g&k J Vj(x)ip 2 k {x)dx 

= ~9k + O (exp (-2\xi - x 2 \/r l )) , \x 1 - x 2 \ » r h 
and its off-diagonal entry is 

(fi,H {2) <f 2 ) = -gi{(P!,(pz) + {<Pi,V2(p 2 ). 
Since v is of finite range, the first term here decays in \x\ —x 2 \ not faster than the second 

— 1/2 

term. Hence, being interested in distances \x\ — x 2 \ that are much bigger than g 12 , we 
can neglect the second term, i.e., we can use as the off-diagonal entry of the matrix the 
quantity —I(x\ — £2), where 

I(xi - x 2 ) = gi&x, <p 2 ) ~ g 2 (<Pi, <f 2 ) (3-14) 

is known as the overlap integral, and in view of (J3.11j) and (|3.13j) we have 

I(x) ~ I Q e~ lxl/ri , \x\ > r h (3.15) 

with 

Io-\E F \. (3.16) 

We obtain that the two lowest eigenvalue of the two-well Hamiltonian can be found 
as the eigenvalues of the matrix 

-9i - x 2 ) \ , 3 ir) 
-I{xi - x 2 ) -g 2 . J ' 

Assuming that g\ > g 2 > 0, we obtain that the eigenvalues of this matrix are 

= -g-(-l) k -W5 2 + P, k — 1,2, (3.18) 

2 In the appendix, we compute exactly the negative spectrum of for v(x) — —S(x) in the 1- 
dimensional case. The results for the conductivity, obtained from this spectrum, coincide with those 
found by using this approximation. 
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where 

8= £L + », = »LZ» (3.19) 

and the corresponding eigenf unctions of the projection of are 

^i(ar) = (fx(x) cos 9 + (p 2 (x) sin0 , g 2Q , 
ip 2 (x) = —(p^x) sin8 + (p 2 (x) cos9, 

where 

tan# = (3.21) 

We are going to use these formulas in the r.h.s. of (|3.9|1 . keeping there only terms with 
m,n = 1,2, i.e., in fact, the term, corresponding to m — 1, n = 2. It is easy to see 
that the equalities £ F = e{ 2) , £ f + z/ = E {2) imply, in view of (EH5|) - (JB~TH|) . that 

v = 2yJ P(y) + 5 2 , y = x 2 — x\. Hence, by (I3.15|) - ()3.16J) and by the condition v <C \Ep\, 
the values of y, contributing to ([3.9)1 . are bounded below by 

r ( I/ )=r,log^ J (3.22) 

and the values of \8\ do not exceed u/2. Under these conditions the coordinate matrix 
element xf 2 ] in (Q: 

x{ 2) = (a* - x 2 ) -1 + (p- 1/2 - ^ 1/2 ) i / xcp 2 (x)dx (3.23) 

2 a/5 2 + P 2VF+P J 

+ , = / xip x {x)^{x)dx (3.24) 



V^P 

between states (|3.2U)) can be replaced by 

x£> ~ (x x - x 2 ) -1 (3.25) 

2a/<5 2 + J 2 

Indeed, the second term in f!3.23|) can be omitted because its ratio to the first term is of 
the order \og2\Ep\/ <C 1. Besides, the term is zero if ip is even. The relative 

order of the third term is the same as the second one. 

In view of the above we obtain that the two-well contribution (|3.9|) to the a.c. con- 
ductivity is 

<t< 2 \v,E f ) = u [ If^A Ep + u- E {2) )8{E F - E? ) ) l i{ Sx )n{ S2 )dg x dg 2 dy. (3.26) 

We integrate first the product of two 5-functions with respect to g\ and g 2 , taking into 
account that \g x — g 2 \ < v <C \E F \ ~ g x>2 . This allows us to replace fi(gi) and pig?) by 
fi(-Ep), to set 

S=\y/" 2 -U 2 (y), (3.27) 

and to obtain in view of ()2.18j) 

<t (2) (v,E f ) = up 2 (E F ) [ \y\ 2p (y) dy (3 28) 

J2\I{y)\>v ^v 2 -AP{y) 
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Note that the restriction 2|/(y)| > v of the domain of integration in (J3.28|) is because of 
the presence of the two (^-functions in (I3.26|) . i.e., in fact, because of energy conservation. 

In view of the inequalities < v << E F and formulas (j3.15|) - (j3.16|) . we can replace 
the condition 2\I(y) \ > v by the condition \y\ > r(u), where r(u) is defined in (I3.22|) . 

The integrand in ()3.28|) is divergent at the lower limit \y\ = r(z/) and decays exponen- 
tially fast at infinity with the rate 2/ri in view of (J3.15)) . Thus the main contribution to 
the integral is due to a /"/-neighborhood of the lower integration limit. This leads to the 
asymptotic expression 

* {2) {v, E F ) = ^^ rf +2 (log ^ ^ . (3.29) 

where Sd is the area of the d- dimensional sphere. Taking into account relations (j3.16|) . 
and (|3.13|) . we obtain finally that 

a^(u,E F ) = p2 P 2 ^ S « \E F \-^ (log^)^ 1 . (3.30) 



In particular, we have for d — 1: 



vV(E F )^ . 3/2 f,2\E 



a(u,E F ) = ^^\E F \-^ log^ . (3.31 
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These are our versions of the Mott formula (jl.3|) . They will be discussed in more details 
in Section 5. 



4 Correlation Functions 

4.1 Generalities 

In this section we study the following two-point correlation functions: 

Ci(x - y; z/, E) = (p E (x, y)p E+u (y, x)), (4.1) 

and 

C 2 (x - y; u, E) = (p E (x, x)p E+v (y, y)), (4.2) 

In writing the above expressions, we took into account the translation invariance in coor- 
dinates of the correlation functions, following from the translation invariance in mean, a 
fundamental property of disordered systems. 

The function C\ of ()4.1|) is closely related to the a.c. conductivity. Indeed, recall the 
spectral theorem, according to which 

Pe(x,v) = J 5{E- E')^ E ,(x)^ EI (y)dE>, (4.3) 

where the symbol J ...dE denotes both the integration over the continuous spectrum and 
the summation over the point spectrum. 

Formulas (jOjl . (I3~4l . and (gZED imply that 

a{u,E) = -— / \x\ 2 C l {x,E,u)dx. (4.4) 



13 



The function C 2 of (|4.2j) is the local DOS-DOS correlator and is a characteristic of local- 
ization, providing information on correlations of eigenstates whose energy difference is v 
and that are localized in spatial domains of distance x — y. 

Comparing (jl.lj) for I = 2 and (|4.1j) and (j4.2j) . we obtain the equalities 

u, E) = K 2 (0, x; x, 0;E,E + v), 
C 2 {x; v, E) = K 2 (0, 0; x, x;E,E + v). { ^ 0) 

We list below certain properties of C\ and C 2 . 
(i) 

\C x {x;v,E)\<C 2 (x;v,E). (4.6) 

The inequality follows from the inequality \p E (x,y)\ 2 < p E (x,x)p E (y,y) that is a simple 
consequence of the Schwarz inequality (ab) 2 < (a 2 ) (b 2 ) and of the spectral theorem ()4.3|) . 



C 1 (x;v,E)dx = 6(v)p(E). (4.7) 

This relation follows from (|4.1|) and (jl.2|) and can be interpreted as a weak form of the 
decay of the correlator C\ at infinity. 
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lim |^V| 1 / C 2 (x;u,E)dx = p{E)p{E + v). (U 

A->oo J A 



To prove this formula, we use the ergodic theorem for p E (x,x), implying the validity of 
the relation 



lim \A\-' [ p E (x,x)dx=(p E (0,0))=p(E) 



on almost all realizations of the random potential. Applicability of the ergodic theorem 
follows from the translation invariance in mean and the decay of the spatial correlation 
in disordered systems (see e.g. [2T]). 

Formula ()4.8|) expresses the decay of correlations between two density operators in 
(14. 2|) as \xi — x 2 \ — > 00. Indeed, its r.h.s. is the product of the averages of these two 
operators (see (jl.2j) ). and its l.h.s. is a weak form of the relation lim^oo C 2 (x, E,E + v) = 
p(E)p(E + v). 

(iv) Assume that for a certain E 

C a (x; v, E) = 8{v)p a {x; E), a = 1,2. (4.9) 

Then 

(a) 

pi(x; E) = p 2 (x; E) = p(x; E) > 0; 

(b) 

P^ E ^ = /^^-^(O^V (4.10) 
where the symbol Ylioc denotes the summation over the localized states only; 
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(c) if one defines the density of localized states as 

Pioc(E) = J p(x ] E)dx=(^5(E-E j )^(0)\, 

then 

p loc (E) < p(E), (4.11) 

and the inequality p loc (E) > is equivalent to the existence of localized states in 
a neighborhood of E, and the equality p loc (E) = p(E) is equivalent to complete 
localization in a neighborhood of E. 



The above properties follow from the spectral theorem (|4.3j) . The functions p a (x; E) of 
()4.9j) are the "diagonal parts" of the r.h.s. of equalities ([4. 5 J) , viewed as functions of two 
variables E\ = E and E 2 = E + v. 

The property (iv) will not be used below. We presented this property to demonstrate 
usefulness of the correlators C\ and C 2 in the theory of disordered system. In particular, in 
the classic paper by P. Anderson jl] the positivity of J p(0; E)dE was used as an indicator 
for localization. The quantum mechanical meaning of J p(0; E)dE is the probability for 
a particle to be in an infinitesimal neighborhood of the origin at time t = 00, provided 
that at t — it was at the origin (the return probability density) |21|. 



4.2 Computations. 

To apply the density expansion formula ()2.10|) to the correlation functions ()4.1|) and 
we write them in the form of extensive quantities per unit volume: 

C a (x) = |A| -1 $ a (a:), a = 1,2, (4.12) 

where 

$i(x) = f C\{{x + a) — a)da = [ (p E (a, x + a)p E+v {x + a, a))da, (4-13) 

and 

&2(x) = I C*2((x + a) — a)da = I (p E (a, a)p E+u (x + a, x + a))da, (4-14) 



E+i 

A J A 



Now it is clear that the role of the functions Fi in ()2.10|) for C a will play $ a , written for 
the /-well Hamiltonian ()2.13|) . 

By using these formulas, the zero- well and the one- well contributions to C a , a — 1, 2 
are absent by the same argument as for the conductivity. The two-well contribution C\ 
to Ci is (cf (COty ): 

C[ 2 \x; v, E) = ~ J] J 5(E + v - E m )5(E - E n ) 

m,n 

x ^ m (a)^ m (a + x)il) n (a)ip n (a + x)p(gi)p(g 2 )dadydg 1 dg 2 , (4.15) 
where y is the separation between two wells, implicit in ij) m n and in E m ^ n . 
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Since v > 0, the diagonal part Yl m =n °f the double sum is zero. Moreover, as in the 
case of the conductivity, we restrict ourselves to the two lowest levels of the spectrum of 
of ipnnj) . found in the previous section in the framework of the projection method. 
This leaves the term m = 1 , n = 2 in the double sum of (j4.15|) . By using (j3.18|) , we can 
integrate with respect to g± and gi the product of two 5-functions, fixing g\ and #2 by the 
relations E + g + ^/ 5 2 + 1 2 = 0, E + g + v — \/ 5 2 + 1 2 = . In view of the condition 
< v <C E we obtain, replacing fj,(—gx) and //(— <7i) by p(-E') in view of (j2.18|) (cf ()3.28|) ): 

Cf } (a;;z/,£) = p 2 (^)^ / daip^a^a) (4.16) 



V»i(a + x)4j 2 {a + x) U =dy. 

2\i( y )\>u ^v 2 -AP(y) 



According to the previous section, if v <C E, the restriction 2\I(y)\ > v is equivalent to 
M > r iy\ where the resonant radius r{y) 3> r\ is defined in (I3.22J1 . This and the form 
(|3.20|) of the functions tp l 2 imply the that if \g± — gi\ < v <C g ~ then 

^(0)^2(0) = costfsintf ^ 2 (a) -^ 2 (a + y)] + O ( e - 2rM/ri ) , 

where u = X\ — x 2 . The formula and the analogous formula with a, replaced by a + x, lead 
to the following asymptotic expression for the two-well contribution C\ to the correlator 
Ci: 

Cf(:r;z/,£) = / (4-17) 



x 



Similar arguments show that the two-well contribution 



[f 2 (a + x) — p 2 (a + x — y)] dy. 



Ci 2 \x;is,E) = ^Yl J ' S(E + u-E m )6(E-E n ) 



(4.18) 



'.2 / \„/.2 



x ^m( a )^m( a + x)/j,(g 1 )n(g 2 )dadg 1 dg 2 dy } 
to the correlator C2 is with the same accuracy: 

ci 2) (x;is,E) = c[ 2 \x;is,E) (4.19) 

+ p 2 {E) / ^ 2 (a)da / —=JL===ip 2 {a + x-y)dy. 
J J\y\>r(v) - 41%) 

We formulate now several properties of C*j 2 ^ and C 2 2 \ following from (|4.17jl - (|4.19jl . 
According to (|4.17|) 

C[ 2 \x;u,E)dx = 



This relation is in agreement with the exact sum rule (|4.7jl . because formula (|4.17jl was 
obtained under the assumption that v > 0. 
Likewise, we have the limiting relation 

C {2) ^p 2 {E), x -> 00, (4.20) 

which is in agreement with the exact sum rule ([4.8)1 . 

(2) 

It is also easy to see that C[ (x; i/, £"): 
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(i) has a positive peak of order 



^(^(log^)"" 1 (4.21) 

at the origin; 

(ii) decays exponentially fast with the rate 2/ri for |x| ^> r/, and is exponentially small 
in the spatial domain r; <C \x\ <C r{y) = 77 log2/ /^ S> rj; 

(hi) has a negative peak of the same order of magnitude (|4.21j) at a ^-neighborhood of 

|x| = r(u); 

(iv) decays exponentially for \x\ ^> r(v) with the rate 2/n, thereby complementing (|4.8jl . 
This behavior of allows us to obtain the Mott formula ()3.29|) from the relations 

(jOD and (mm) . 

The correlator has the same behavior as Cf^ in x till |x| < t(^), in particular it 

(2) 

is exponentially small in x if r; <C |x| r(z/). Then C| becomes asymptotically equal 
p 2 (E) in the domain \x — r(u)\ < ri and it is equal to p 2 (E) for all x, \x\ ^> r(i/) (see 
(|4.20Jl ). In view of spectral theorem one can expect p E (x,x) to be proportional to iI> 2 e(x) 
in the strong localization regime (cf (|5.18Jl ). Then the factorization property (|4.20|) can 
be interpreted as the statistical independence of the localized states of energies close to 
each others with separation much bigger than r{y). On the other hand, the exponential 
smallness of for T\ ^> \x\ ^> riy) can be interpreted as a kind of strong correlation 
between states close in energy, that are not sufficiently well separated in space. These 
correlations can be viewed as a manifestation of a certain "repulsion" of nearby levels 
in the sense that the probability that nearby levels are close tends to zero as the level 
spacing tends to zero (see [121 El) for discussions of this property). Figures 1 and 2 show 
examples of graphs of Ci t 2, 




Figure 1: One dimensional correlation function C\{x) and C*2(x) with v = 10 4 , r; = 1 and 
p{E) = 1. 
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Figure 2: Two dimensional correlation function C±(x) and C2(x) with v = 10 4 , 77 = 1 and 
p(£) = 1. 



Results similar to those outlined above were obtained as asymptotically exact ones in 
[T4*j in the one dimensional case of the strong localization regime and in ^QJ in the 
one dimensional case of the weak localization regime (see the following section for more 
details). We see, however, that in dimension greater than 1 the characteristic value (j4.21j) 
of the peak of the correlation function (|3.2|) diverges as the difference of the energies 
tends to zero. Thus, unlike the conductivity that has the log-factor in all dimensions, the 
correlation functions C\p are logarithmically big in the energy difference for |x| ~ r{y) 
only in dimension bigger than 1. Similar results were obtained in jT^j in the frameworks 
of the instanton approach (see Subsection 5.4). 

The "two- hump" states ()3.20|) appear in our approach just as a computational tool, 
allowing us to find leading contributions to the low frequency conductivity and to the 
correlators C\ t 2, by using the density expansion of Section 2.1, just as the "one-hump" 
states ()3.11|) are necessary to find the low energy asymptotic of the density of states in 
our approach (see also Section 2.3), in the optimal fluctuation method [2H1I2I1IZ1, and its 
version, known as the instanton approach (see Section 5.4, and references therein). 
On the other hand, the development of localization theory of the last decades suggests 
that the "one-hump" states carry certain information on the structure of genuine localized 
states in disordered systems. This suggests the belief, that the "two-hump" states also 
reflect certain properties of genuine localized states. If this is true, we can interpret the 
above results on the spatial behavior of the correlators C\^. in the following way. The 
existence of the length scale r[y) of ()3.22|) . that determines drastic changes of the spatial 
behavior of the correlators is due to the "interaction" between close energy levels, 
and the interaction mechanism is the resonant tunnelling between the "bare" one- hump 
states, i.e., between different centers of genuine states. The parameter I of ()3.15|) - ()3.16|) 
is the characteristic interaction energy, determining the level splitting (spacing), and r(z/) 
is the tunnelling distance, determined by the two energy scales (E = (Ei + E 2 )/2, v = 
\E2— E\\, E » u). This inter- level interaction is a mechanism of a certain level repulsion, 
that prevents the spatial domains where the states are essentially non zero to be close 
and, as a result, leads to the exponentially small values of the two-point correlators for 
r\ « \x\ « r(v). 
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5 Discussion 



5.1 Corrections 

We comment now on the corrections (next terms of the density expansions) to our formulas 
of Sections 2-4. We are not able to prove the convergence of these expansions. We simply 
argue that they should be asymptotic, i.e., that their terms should be small in successive 
powers of p(E). We will begin from the density of states itself. 

It is easy to see that the next term in the expansion of the DOS has the form 



P (2) 



(E) = j { [s(E - E[ 2) ) - 5(E 

+ [5(E - E®) - 5(E - 4 2) )] } p{ gi )p{g 2 )dydg x dg 2 . 

where E® are given by and e? = -max(^ 2 ), ef = -min( gi ,g 2 ). Recall that 

we assume that p(g) is smooth enough and decays sufficiently fast for large g. Thus p^ 
will be of the order 0(p 2 ) if the integral in the relative distance y between the wells will be 

convergent. This fact follows from the inequality \E k — e^\ < (y 8 2 + P(y) — \8\) < \I(y)\, 
the exponential decay of I(y) (see (|H.15jl ). and the smoothness of p(g), allowing us to 
transfer derivatives of delta-functions to the /x's. 

In the general case of the correction of the order I the appropriate integrals in relative 
distances between wells will be convergent because of the subtractions of the functions 
Fk of lower orders k < I from that of the order I in the Ith term of the density expansion 
(|2.1()j) . the sufficiently fast splitting (additive clustering) of negative eigenvalues E^ of 
the /-wells problem into the sums of negative eigenvalues E^ k > of the fc-wells problems 
k < I and again because of the smoothness of p(g). 

The situation is less simple in the case of the conductivity as we have seen already 
for I = 2. This is because of the presence of families of tunnelling configurations for any 
number of wells (for example, for / = 3 there are two families: the equilateral triangles 
and the three equidistant points on a straight line). These configurations are responsible 
for the absence of decay (and even for the polynomial growth) in distances between 
wells of matrix elements xfj on the corresponding resonant sub-manifolds and for the 
appearance of extra powers of logz/o/V (where vq can be different from that of formula 
()3.30p ). However, since the dimension of these resonant manifolds grows slower than /, 
these powers of logu^/v will be always multiplied by powers of v, given by the dimensions 
of the manifolds transversal to the resonant ones. This is why the higher terms in the 
expansion of the low frequency conductivity should be small compared to the terms in 
Mott's formula (PI) . 

In other words, it seems reasonable to believe that these higher resonant configura- 
tions will produce new peaks and new length scales in the higher terms of the density 
expansion of the correlators, but that the amplitudes of the peaks will be small relative 
to the amplitude (|4.21|) of the peak due to the resonant pairs. One can also speculate 
that for bigger densities of states (i.e., for energies closer to the mobility edge) higher 
resonant configurations will play a more significant role, leading eventually to the loss of 
the exponential decay of the correlators and to the derealization transition according to 
the scenario, outlined in EOj 

5.2. Asymptotically exact one- dimensional results 
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The asymptotic behavior of the low frequency conductivity in the strong localization 
regime of the one- dimensional Gaussian white noise potential, defined by the relations 



(V(x)) = 0, (V(x)V(y)) = 2D5(x - y), (5.1) 

was studied in [T3j. The potential is often used in the theory of one- dimensional disordered 
systems (see for results and references). In particular, the density of states p(E) and 
the Lyapunov exponent ~f(E) of the Schrodinger equation with this potential can be found 
in quadratures. The strong localization regime corresponds to negative energies of large 
absolute value 

L> 2/3 <|£|. (5.2) 
In this case we have the following asymptotic formulas [21] 

p (E) = &e-^ /5 / 3D , >y(E) = |£| 1/2 - (5-3) 
nD 

Moreover, the rate of the exponential decay of the eigenf unctions ip E is j(E), because we 
have with probability 1 j^U HZj : 

lim Ixl-Mog^Kx) +^|(x)) 1/2 = - 7 (£) (5.4) 

\x\— >oc 

Hence, the exact asymptotic form of (|5.3|) for the localization radius 

n(E) = l/ 7 (£). (5.5) 

coincides with our approximate formula (|3.13|) . 

In the paper (Hj the low frequency conductivity was found using the Grassmann 
functional integral representation of the Green's function, which leads to an integral rep- 
resentation for the correlator C\ of (|4.1|) (recall that the conductivity is related to the 
correlator via formula (|4.4jl ). The condition ()5.2j) allowed the authors to apply the saddle 
point method to this integral representation. We will summarize the results of J3] in a 
form close to that of Sections 3 and 4. 

The "two-hump" states, similar to (|3.2(J|) appear in [14 as the saddle points of the 
effective action for C\. The states have in general a rather complicated (two-instanton) 
form, but in the low frequency limit < v <C \E\ they can be written in the form (|3.20p . 
in which the role of the "bare" states play 

< ^ 1 ' 2 ^ V2r~i cosh (x ± y/2) jri ' 

where y > y^v) and 

yv{y)=r l \og%\E\/v (5.7) 

(cf (J3TTJ), (jH32D, and As for the angle 9 of (jHHOD , it is defined by the relation 

e -y/ri = e -yo/r t g j n th&t can be written as 

e -y/ri 

tan^ = . (5.8) 

g-2/O/n -)_ ^Q-tyo/rt _ e -2y/r l 
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Introduce I(y) = I$e y ^ n , where Jo = 4|-E|. Then formula (J5.7)) can be written as I(yo) = 
u/2. These formulas have to be compared with ()3.15|) . and (j3.16|) . Furthermore, setting 



5 = yJuyA-P(y) = yJp(y )-P(y), (5.9) 

(cf (|3.27|l ). we can write (|5.8J) in a form, analogous to that of (j3.21j) . 

According to [13], the correlator C\ has the following asymptotically exact form for 
< v < \E\: 

/f e~ y l Tl 
da / + x)il) 2 (a)ip 2 (a + x)—= =dy, 

Jy> yo ' ' y/e-zyo/n _ e -2y/n 

which can be written as (j4.16|) because, in view of the above notations, we can write the 

expression e~ y l Tl (e~~ 2y °l ri — e~ 2y ^ Vl ) ^ 2 in the last formula as v {v 2 — AI 2 (y)^j 

Likewise, the asymptotically exact expression for the low frequency conductivity, ob- 
tained in [Tl] . coincides with our formula f)3.28|) . and the correlator C 2 has the form (|4.19j) . 
after the replacement Ep — > AEp under the log sign. The correlator C 2 was not considered 
in jHj, however it can be found by using the techniques, developed in the paper. 

We note a certain difference of these asymptotically exact results and our results. 
Namely, the role of the resonant distance r(u) of ()3.22|) in the results of J3] plays (15. 7|) 
that differs from ()3.22|) by the factor 4 under the logarithm. A possible simple reason 
for this difference can be the fact that our estimate (j3.16|) for the amplitude Iq of the 
overlap integral indicates only its order of magnitude, but not its precise value, or, more 
generally, that the projection method is not precise enough. 

5.3. Weak localization regime in one dimension. 

The case of the Gaussian white noise (J5.1|) in one dimension has also been studied in the 
weak localization regime of large positive energies 

D 2/3 < E (5.10) 

(see the works (HI d EH [TTJ, EH]). The density of states in this case is the free one 
p (E) = (2tcE)~ 1 ^ 2 , and the localization radius is 

AE 

n = -^- (5.11) 

The rate of the exponential decay of wave functions is 1/rj with probability 1, as it was 
in the strong localization regime (see (|5.4jl ). 

There are several techniques that can be used in this case d ED EL an d yield 
the low frequency conductivity and the correlators C\ and C 2 in quadratures. It turns 
out that these quantities have qualitatively the same spatial behavior as in the strong 
localization regime, provided that 2E F (i.e., 2I according to (|3.16Jl ) in (|3.31|) is replaced 
by (D/2E]I 2 ). Note that (D/2E l J 2 )~ l coincides with the relaxation time r, well known 
from the kinetic theory [T0| . According to JS], the quantities Iq and r -1 have the same 
meaning: they give the order of magnitude of the difference of the energies (spacing) of two 
localized states, whose centers are separated by a distance of the order of the localization 
radius. Similarly, the role of the resonant distance in the two-point correlators C\ )2 plays 
(cf (1T221 and (JST}): 

r{y) = r/log8/z/T, (5-12) 
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and the rate of the exponential decay of the two-point correlators G\i near the origin 
is l/2r/. This rate is 4 times less than the rate 2/77 of these correlators in the strong 
localization regime, found in Section 3 from the naive prediction, based on the spatial 
behavior of the envelope of the eigenf unctions with probability 1 (see (|5.4jl ). and in [T4~j 
from an asymptotically exact analysis of the corresponding correlators. This difference 
can be related to the fact that eigenfunctions in the one dimensional case in the weak 
localization regime are much more spread out than in the strong localization regime. 
Hence, their behavior on almost all realizations can differ from the behavior of their 
moments, entering exact formulas ([3.2)1 . ([4.1)1 . and ([4.2)1 . 

We stress that the basic properties of the strong localization regime and, in particular, 
those, motivated assumptions and techniques of this paper, are different in several impor- 
tant points from the basic properties of the weak localization regime in dimension 1, where 
the mechanism of localization is not trapping in deep and rare localization wells but the 
enhanced backscattering due to the destructive interference between incident and reflected 
waves from many defects. One of manifestations of this complex statistical structure of 
wave functions in the weak localization regime is the value of the rate of exponential decay 
of the correlators Ci^, discussed above. Moreover, according to [TT], the characteristic 
length scale of the correlators C% t 2 in the neighborhood of r(v) is y/HjTjri, i.e., is much 
bigger than the scale r\ in the neighborhood of the origin, while, according to our formulas 
and respective formulas of jT3] , in the strong localization regime this scale is r\ both near 
r[y) and the origin. 

5-4- Instanton approach. 

This is a version of the variational method, proposed first by I. Lifshitz j2U] to find the 
asymptotic form of the density of states and other characteristics of disordered systems 
in the strong localization regime (see e.g. [7j). The instanton approach was used to 
analyze the correlators Ci^ and the low frequency conductivity for the white noise ran- 
dom potential in d dimensions in paper ^5], in which the reader can find references on 
earlier applications of the approach. It is based on the assumption that in the strong 
localization regime the two-point correlators correspond to the two-well potential that 
minimizes the total probability distribution of the random potential under the constraints 
that H{V)ip k = Ekip k , k = 1, 2 and that the well centers of the "optimal" potential are a 
distance y = x\ — x<i apart. This has to be compared with the DOS computation, where it 
is assumed that the optimal potential is a well for which H(V)ip = Eip (see |21) I7j). The 
derivation of final formulas in [Tj£. is rather involved because of existence of two energy 
scales and of collective modes, in particular those that correspond to the center of mass 
(xi + X2V2 of the optimal potential (it is an analogue of our parameter a in (j4.12[) - 
()4.14)l ). As a result, it is shown in ^5] that in the strong localization regime (called the 
hydrodynamic regime in [T^]) the correlator C\ and the low frequency conductivity have 
qualitatively the same form as those found in Sections 3-4. 

We note also that the 1-dimensional results for the white noise potential of ^3] can be 
viewed as a justification of the instanton approach in the one-dimensional case, because 
it was shown in this paper that the two-well potential of a special from is indeed a saddle 
point of the respective functional integral. 

5.5. Maryland model. 

The most widely known signature of localization is the exponential decay of the localized 
states at infinity. However, the initial derivation of the Mott formula ([1.3)1 as well as 
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the above derivation are based not only on the exponential localization, reflected in the 
exponential decay of "bare" states of the independent quantization in each localization 
well, but also on the weak correlation between the spectra of independent quantization, 
reflected in statistical independence of localization wells in our effective potential (J2.1J) 
and in appearance of the "two-hump" states in our calculations of Sections 3, and 4. 
The relevance of the last property becomes clearer if one recalls the results, obtained for 
an explicitly soluble model of an incommensurate system, known as the Maryland model 
^3JI2nil2E]- This is a multi-dimensional tight binding model with an arbitrary short-range 
and translation invariant hopping and with the potential of the form 

V(x) = g tan 7r(a ■ x + u), x & Z d , (5.13) 

where g > is the coupling constant, a is a (^-dimensional vector with incommensurate 
components, and u G [0, 1) is a phase, that plays the role of a randomizing parameter. It 
was found in the mentioned papers that if for some C > and (3 > d the vector a satisfies 
the Diophantine condition 

\a-x + m\ > C/\xf (5.14) 

for any integer m and i^O, then all the states of the model are exponentially localized 
for any coupling constant, any energy, and arbitrary dimensionality d of the lattice Z d . 
Since the potential has arbitrary high peaks, the model can be viewed as an explicitly 
soluble model of the strong localization regime. The spectrum of the model consists of 
the solutions of the equation 

N(E t (u)) = a -t + uj (mod 1), (5.15) 
where t is a lattice point, N{E) = p(E')dE', 

p{E) = - [ 9 —„ dk, 

vr J Td (w{k) - E) + g 2 

is the density of states, in which w(k) is the Fourier transform of the hopping coefficient, 
and T d is the (i-dimensional torus. 

It is easy to show that for each point t of the ci-dimensional lattice the equation has 
a unique solution, that if E tl {uj) = E t2 (u), then ti = t 2 , and that the set {E t (uj)} t&Z d of 
eigenvalues is dense for any uj e [0, 1). 

The corresponding eigenfunctions if> t , t G Z d have the form 

il> t (x)= X (x-t,Et(u)), (5.16) 

where x{ x i E) decays exponentially in x: 

\x(x,E)\ < Ce~ lxl/Tl (5.17) 

with some positive ri(E). Formulas ()5.15|) - (j5.17|) seem fairly natural in the case of the 
strongly incommensurate potential (|5.13|h where due to the absence of any symmetry the 
only good quantum number to label levels and states is the "center" of localization well. 

One can also say that Mott's notion of the localization centers is explicit here, because, 
according to (|5.15j) and (|5.16|) . for any lattice point t there exists a unique eigenvalue 
E t , whose eigenfunction is exponentially localized in a neighborhood of t. Thus the set 
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of localization centers coincides with the whole lattice and the density of localization 
centers, whose states have energies in a neighborhood of a given E is the density of states 
p{E). This fact can be interpreted as the uniform distribution in space of the localization 
centers, corresponding to energy E, and is in qualitative agreement with our assumptions 
of Section 2, formula (J2.18j) in particular. 

On the other hand, the low frequency conductivity and the correlators C\ and C2 for 
the potential ()5.13|) have a rather different structure than in the case of random potential 
discussed in Sections 3 and 4. This can be seen from the form of the kernel p E (x,y), 
following from (j5.15|) - (|5.16j) : 

Pe(x, y) = J2 5{ * E " E ^)Mx - *, E t (u))x(y - t, E t (u)). (5.18) 
tez d 

Consider first the correlator G%. Plugging (j5.18|) into ()4.2|) . and recalling that the averaging 
operation (...) here is the integration with respect to the parameter to G [0, 1) of (|5.13jh we 
find first of all that the correlator C2 is not a regular function. Rather, there exists a dense 
set of special frequencies for which C2 has 5-peaks. If, however, we are interested in the 
gross features of C 2 , then we can apply a certain smoothing procedure, say z/ _1 ...dp' . 
Then we obtain that there exists the length scale 

nM = {^p-) > ME) = C/p(E), (5.19) 

(here (3 and C are defined in (|5.14Jl ). such that C 2 (x; u, E) is of the order e~ ri( - u ^ ri if 
\x\ <C and C2{x;u,E) is p 2 {E) if \x\ » T\{y), and the transition from the first 

value to the second one is in the layer \x — ri(v)\ ±2 rj, where r\ is defined in (|5.17j) . 
We see that the qualitative form of the correlator C2 for \x\ >> r\ is similar to that in 
the random case, however there is no peak at the origin and the length scale ()5.19j) is 
polynomial in v (cf (|3.22|) ). In addition, the length scale ()5.19|) has a different origin than 
(13.22)1 : it is not due to the tunnelling for "soft" resonant pairs, but due to the Diophantine 
condition ()5.14j) . which determines now the distance to the nearest localization well of an 
almost same energy. At low frequencies r\(v) is much bigger than the resonance tunnelling 
distance r{v) of ()3.22j) . This leads to the qualitative change of the form of the correlator 
C\. Indeed, by using the same argument, we find that G x is of the order e~ 2r ^/ ri < 1 
for all x. This and formula (|4.4)l imply that the low frequency conductivity is of a similar 
order 

a(u, E F ) - exp{-(z/ 1 ( J E;)/z/) 1 /' 3 }, v x = (5.20) 

r i 

The significant difference between (I5.20|) and (jl. 3)) can be related to the absence of long 
range tunnelling in the Maryland model. The spectrum of the model is too "rigid", the 
energy levels are too regularly distributed and small level spacing are too rare for the long- 
range tunnelling to happen. This illustrate the role of resonance tunnelling in obtaining 
the Mott formula as well as the range of applicability of the approach of this paper, based 
on ansatz (|2.1|) and on the density expansion. Besides, we see that the low frequency 
conductivity provides a physical distinction between the strong localization regimes of a 
random shortly correlated and smoothly distributed potential, and the incommensurate 
potential ()5.13|) (recall that the density of states and the Lyapunov exponent coincide 
for the Maryland model and for the random model in which the potential is a collection 
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of independent identically distributed Cauchy random variables, and in which we expect 
our approach to be applicable). Besides, recalling the structure of the localized states 
for smooth incommensurate potentials of large amplitude in one dimension [2*5] . e.g. the 
potential g>cos27r(ax + o;), g » 1, one may expect that these potentials will be closer to 
random potentials in the spatial behavior of two-point correlators and the low frequency 
asymptotics of the conductivity. 



Appendix 

One-dimensional case with delta potentials 

To support the usage of the projection method by which the bound states of the two-well 
Hamiltonian in Section 3.2 were found, we will consider here the one-dimensional case 
with two delta-wells. The corresponding Hamiltonian is: 

H {2) = ~ - 2y/gi6(x - x x ) - 2^g- 2 5(x - x 2 ), (A.l) 

where gi }2 > 0. In this case each of two one-well Hamiltonians 

d 2 



2^MJ5(x - x 12 ) 



dx 2 

has the unique bound state 

Vi,2( x ) = 9\[iv(y/9i? x )i = e ~ W > ( A - 2 ) 

corresponding to the energy 

E§ = -g lt2 . (A.3) 

Since the Hamiltonian is invariant under translation, we can replace X\ by 0, and x 2 
by y. It is easy to see that has two bound states: 



^1,2^) = vW(o) ex P (-VW\\ x \) + sfgi^iv) ex P \-VW\\ x - y\) » ( A - 4 ) 

\ / \ /J E=E\ 2 

(2) 

where E\ ^ are the corresponding energies. They solve the equation: 



(vW\ - Vg~i)(V\E\ -V&) = v^ ex P (-^VWWyl) (A.5) 



Assuming the same accuracy as in Section 3.2 (\gi — g 2 \ g± t 2, \y\ » 91,2)1 we find 
that the solutions e[ 2) 2 of |A~5|l have the form (j3~TI"jl in which I(y) = 2ge~^ yl (cf flTToT) - 
f|3.16)l ). and the eigenfunctions (|A.2|) have the form ()3.20j) - f|3.21j) in which <f 1 2 are given 
by (CH). 

Another way to act in this case is to plug the exact states and levels, given by (|A.4|) - 
(IA.5J) . into the expressions ()3.9jl . (|4.15jl and 1)4.18)1 for the two well contributions for the 
conductivity and the correlators C\, and C 2 . This leads to rather complicated formulas 
which, however, have the same asymptotic behavior as our formulas (j3.28j) . (j4.17|) . and 
f!4.19|) in the asymptotic regime < v <C \E\. 
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